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Using the example of the SU (2) gauge theory in 3 + 1 dimensions we consider the construction of 
a 3-dimensional effective model in terms of Polyakov loops. We demonstrate the application of an 
equilibrium self-consistency condition to the systematic analysis of the contribution of various (global 
Z(2) symmetric) terms in the effective model action. We apply this analysis to the construction of 
a simple effective action with the minimum necessary number of operators. Such an action is shown 
to be capable of reproducing relevant observables, e.g. the Polyakov loop ensemble average, within 
the desired accuracy. 
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I. INTRODUCTION 

There has been renewed interest in the construction of 
effective models for the finite temperature confinement- 
deconfinement phase transition of non-abelian SU(N) 
gauge theories, see Ref. [l[ for SU(3) and Refs. @,|| ror 
SU(2). The idea is to use relevant variables to rewrite the 
gauge theory as an effective lower dimensional spin model 
with the goal of simplifying the description of phase tran- 
sition dynamics. For SU(N) gauge theories the center 
group global symmetry Z(N) controls the transition: it 
is spontaneously broken in the deconfincd phase but un- 
broken in the confined phase. Therefore, the effective 
model can be expressed in terms of the relevant order 
parameter (Polyakov loop) [J, Q : 

Nt j 

^x = n^+nt> L(x) = -TrV x . (1) 
n=l 

For SU(2) gauge theory the transition is second order 
and the universality argument holds, ft is conjectured 
that for N > 2 groups the integration of all degrees of 
freedom but the order parameter leads to a short-range 
effective model 

The successful construction of such an effective model 
results in significant simplification of the description of 
the physics. Also, in practical terms, it is much easier to 
perform Monte Carlo simulations of a (d — l)-dimensional 
spin system than a d-dimensional gauge theory. 

Using pure symmetry considerations it is straightfor- 
ward to identify an infinite set of possible operators that 
may enter in the effective model construction. As the 
next step, it is necessary to order these operators by their 
importance. This would allow for a discriminative trun- 
cation of the effective action to a manageable minimal 
set of operators. It is possible to use the strong coupling 
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and character expansions, together with natural trunca- 
tion criteria, to sort the hierarchy of operators of a given 
interaction range 0j- The interplay, however, of truncat- 
ing in the coupling versus the range of hopping terms in 
the action cannot be resolved by this approach. 

In general, truncations result in an action which is ca- 
pable only of an approximate description of the orig- 
inal gauge theory ensemble. Passing to the effective 
model defined by the truncated action may introduce 
non-equilibrium effects (thermalization) . This is espe- 
cially important since the effective theory parameters 
are estimated starting from the gauge theory ensemble, 
which may not be a representative equilibrium ensem- 
ble of the effective action. Therefore, it is important to 
control residual non-equilibrium effects. 

The most interesting aspect of the study of non- 
equilibrium effects, however, is that it provides a natural 
classification of the importance of various effective model 
operators. This is achieved by adhering to the following 
procedure: 

• From a given ensemble of 4D gauge configurations 
U x .x produced by standard gauge theory methods 
(heatbath/overrelaxation), we extract a 3D field of 
Polyakov loops V x . 

• An effective action is chosen by specifying a partic- 
ular set of operators. 

• The system is allowed to evolve microcanonically 
along the sheet of constant energy measured ac- 
cording to the effective action. 

• The process of equilibration is observed. Particu- 
larly we measure the evolution of various physical 
observables as they thermalize at new values. 

The microcanonical evolution is realized through de- 
mon updating 8, 9]. We conduct our study over a range 
of lattice couplings. Since, however, the normalization 
procedure for the demon energy is unknown we decided 
to study only the ensemble average of the fundamental 
Polyakov loop (L). Its evolution is monitored as the 
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system is evolved toward equilibrium. We quantify the 
overall non-equilibrium effect by the difference between 
the final (equilibrium) value and a starting (gauge en- 
semble) value 6(L). If one could chose an effective ac- 
tion capable of reproducing the dynamics of Polyakov 
loops exactly, then no discernible change in observalbes 
would be observed. By adding or removing operators 
and then observing the change, we can sort these opera- 
tors by their importance. The truncation criteria can be 
set so as to minimize the non-equilibrium effects (ther- 
malization flow) in the transition to the effective model, 
starting from an ensemble generated with the gauge the- 
ory Lagrangian. This approach is similar to equilibrium 
self-consistency ideas used in previous MCRG decimation 
studies OH Hud- 

Note that the use of the demon method allows for 
easy monitoring of observables while performing the mea- 
suremts. This is in contrast to previous studies rely- 
ing on the Schwinger-Dyson method for measurements 
of effective couplings. In that method "black box" mea- 
surements are performed without any control of non- 
equilibrium effects. One could improve on the SD proce- 
dure by performing the Monte-Carlo evolution with the 
new measured action starting from the gauge configura- 
tions, while observing the equilibration of observables. 
Then the analysis of effective action terms similar to the 
one proposed here should be used with the end goal of 
avoiding any equilibration changes in observalbes. 

It should also be noted that the equilibration effect can 
be easily missed on small lattices: 20 3 x 4 shows almost 
no effect. 

This approach will be pursued in this paper in the 
case of the 577(2) gauge theory. The same method is 
also applicable to other SU(N) gauge theories, and, in 
particular, SU(3) gauge theory. 



II. THE SIMPLEST ONE-OPERATOR 
EFFECTIVE ACTION 

It is natural to start building an effective model consid- 
ering only the simplest possible term. Therefore, in this 
section we focus on the fundamental character nearest 
neighbor hopping term 
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After preparing a typical gauge configuration at gauge 
coupling (3 we apply the demon method assuming the 
effective action Si = \\A\. The demon was originally 
thermalized for 100 sweeps on an equivalent gauge con- 
figuration which was then discarded. In this way we re- 
move the effect of demon thermalization. In Fig. [T] we 
present the evolution of the Polyakov loop average (L) 
at fixed physical temperature T — 4T C and different lat- 
tice couplings. Since the renormalization of the Polyakov 
loop is multiplicative, in order to compare the relative 
flow scales at different j3 we normilize the Polyakov loop 
to unity. One notable feature of the plot is that the 
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FIG. 1: The renormalized Polyakov loop average at T - 
on different 40 3 x N t lattices. 
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TABLE I: The relative change of the Polyakov loop ensemble 
average 8L/L = 5(L)/{L) after 1000 Monte-Carlo sweeps at 
various temperatures and inverse lattice couplings on 40 3 x N t 
lattices. The starred (*) value is computed using the weak 
coupling 1 loop formula. 





T = 4T C 


T = 2T C 


T = T C 





N t SL/L 


N t SL/L 


N 


3.00* 


8 0.974(4) 




32 


2.7310 


4 0.191(1) 


8 0.973(3) 


16 


2.6355 


3 0.0501(4) 


6 0.975(3) 


12 


2.5104 


2 0.0073(1) 


4 0.57(1) 


8 


2.2991 




2 0.0153(3) 


4 


1.187348 
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Polyakov loop shows vanishingly small change as (3 val- 
ues become smaller. 

In all measurements presented here there are several 
independent simulations which allow us to extract errors 
using the jack-knife method. For this simulation we take 
[3 values to correspond to critical coupling values at some 
integer Nt (see the last column in Tab. [J), which were 
estimated in [l3j . This makes fixing the temperature very 
straightforward. 

In Tab. |I]we present the data for the relative change in 
the Polyakov loop 5(L)/(L) after 1000 sweeps evolution. 
It shows insignificant flow as one approaches the strong 
coupling region while the flow in the the weak coupling 
region is considerable. 

In the limits of strong and weak coupling the consid- 
ered action is expected to reproduce the physics of the 
underlying gauge theory [l4j]. It is possible to estimate 
the coupling of the spin model Ai using perturbative ex- 
pansion. In the strong coupling limit it is [3], 0HH 

Ax = (f3/4) Nt - (3) 
In the weak coupling one expects 

Ai - /3/(N 2 N t ) = (3/{AN t ). (4) 



3 



This part of our study is similar to [15| . We measure 
the effective coupling using the demon method. For this 
we allow the demon to thermalize the system for 500 
sweeps, and then monitor the demon energy value for the 
next 10 sweeps. We present the values of coupling in Tab. 
[UJfor various (3 values and 3 lattices. In Fig. [2]we plot the 



TABLE II: Effective coupling for the simple one-operator 
model as a function of the gauge coupling for N t — 2, 3 and 4 
lattices. Last rows show the leading coefficients of the fits. 






N t = 2 


N t =3 


N t =4 


0.5 


0.0174(4) 


0.0022(2) 


0.0006(2) 


0.8 


0.0430(2) 


0.0082(2) 


0.0016(4) 


0.9 


0.0542(2) 


0.0114(4) 


0.0026(2) 


1.187348 


0.0924(1) 


0.0267(4) 


0.0078(2) 


1.6 


0.1684(2) 


0.0672(2) 


0.0262(2) 


1.8 


0.2236(4) 


0.1008(2) 


0.0456(4) 


2.0 


0.3150(2) 


0.1502(2) 


0.0808(2) 


2.2991 


0.3540(1) 


0.3042(2) 


0.1978(14) 


2.5104 


0.3783(1) 


0.3170(2) 


0.3001(1) 


2.6355 


0.3921(1) 


0.3248(2) 


0.3042(2) 


2.7310 


0.4028(1) 


0.3300(4) 


0.3072(2) 


3.0 


0.4326(2) 


0.3464(2) 


0.3166(2) 


a s 


0.0642(4) 


0.0159(3) 


0.0038(1) 




0.0625 


0.0156 


0.0039 


Q>w 


0.1109(1) 


0.0599(5) 


0.0334(5) 


1/4JV* 


0.125 


0.083 


0.625 



effective coupling values. Note that the figure contains 
more data points than the table. We check the weak 




FIG. 2: Effective spin model coupling as function of gauge 
coupling. Arrows indicate the critical coupling values for lat- 
tices N t = 2 and 4. 

and strong coupling behavior against the perturbation 
theory results. For the strong coupling we perform the 
fit to a s (3 Nt +c, while for the weak coupling the functional 



form is a w (3 + c' . The results of the fits are presented in 
the lower part of Tab. [Til We indicate by solid lines in 
Fig. [2] the perturbation theory results. 

We observe that at strong coupling the perturbation 
scaling is indeed achieved. This nicely corresponds to 
the fact that there is minimal flow in this region. For 
weak coupling we observe deviation from the perturba- 
tive result (similar observations were made in [l5| and 
[3j). This is easy to understand since in this region we 
observed significant flow of the Polyakov loop during the 
thermalization. Note that smaller time extent lattices 
show results closer to the perturbation theory behavior 
in the weak coupling regime. The smallest N t = 2 lattice 
indeed shows fairly good agreement in this region. 

We conclude that the one-operator action is capable of 
describing the gauge theory in the strong coupling limit. 
In the weak coupling limit, however, it represents a good 
approximation only for small time extent lattices N t — 1 
and 2. This result is not surprising since the gauge theory 
at weak coupling exhibits non-local nature which cannot 
be captured by a local effective action. 



III. ANALYSIS OF GENERAL EFFECTIVE 
ACTION 

Next we fix the coupling value to (3 = 2.7310. This 
is a value from the region where, as we saw in the pre- 
vious section, the simplest action exhibited a significant 
flow in the Polyakov loop average. On 40 4 x 4 lattice it 
corresponds to fixed temperature 4T C . We start adding 
higher order terms to the simplest action in order to re- 
duce the flow. It is natural to group the terms by the 
spatial range of the hopping terms. The nearest neigh- 
bor hopping term group y = x + /t is 



il = 


Xl/2(Px)Xl/2(Py), 


0(L 2 ) 


(5) 




Xl{V x )Xl{Vy), 


0(L 4 ) 


(6) 


h = 


Xa/2{V x )xi/2{V v ), 


0(L 4 ) 


(7) 


A 4 = 




0(i 6 ) 


(8) 


A 5 = 


X2(V x )Xl(Vy) 1 


0(L 6 ). 


(9) 



The potential terms group is 



Bi = 


XlPx), 


0(L 2 ) 


B 2 = 


X2(Px), 


0(L 4 ) 


B 3 = 


Xs{V x ), 


0(i 6 ) 



(10) 

(11) 

(12) 



Here the order of operators 0(L n ) is specified by repre- 
senting characters as nth order polynomials of the fun- 
damental representation loop L. The next to nearest 
neighbor group is 



Cl = Xl/2(V x )Xl/2(Vy), 

C 2 = Xi{V x )xi{V v ), 

Ca = X3/2CPx)Xl/2(Py)- 



(13) 
(14) 
(15) 
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FIG. 3: The Polyakov loop flow for the simplest one-operator 
action, various two-operator actions and the A12 + B\ + C12 
action on a T = 4T C f3 = 2.7310 lattice. 

By analogy we introduce Di(x,y) for the next to next 
to nearest neighbor terms y = x + fi + + k and [i ^ 
v 7^ k. We consider various actions of the form Si — 
^OiOi, where O is any of A,B,C or D. For simplicity 
we introduce the shorthand notation Oi...j...k = S\ + ... + 
Sj + ... + Sk, thus A12 = XaxAi + \a 2 A 2 . 

TABLE III: The change of the Polyakov loop L after 1000 
sweeps for various combinations of terms in the effective ac- 
tion. The original gauge theory configuration is at j3 — 2.7310 
on a 40 4 x 4 lattice. 



s 


5L 


A 1 


0.1556(1) 


Ai + Bi 


0.047(1) 


A12 + Si 


0.031(1) 


A13 + Bi 


0.045(1) 


Ai + B 12 


0.044(1) 


A123 + B12 


0.034(1) 


^4l234 + B\2 


0.0307(5) 


A1235 + B\2 


0.0317(5) 


A123 + Bl23 


0.0309(2) 


^12345 + -B123 


0.0286(4) 


Ai + Bi + Cl 


0.0142(4) 


A12 + Bi + Ci 


0.0102(3) 


A12 + B\ + C'12 


0.0065(2) 


Ai2345 + -B123 + Cl23 


0.0041(1) 


Ai + Bi + Ci + Di 


0.0102(1) 


A u + Bi + C12 + Di 


0.0048(1) 


A 12 + Bi + C12 + D 12 


0.0042(1) 



A12345 + B123 + C123 + D12 0.0049(1 



It is impractical to consider all combinations of the 
indicated operators in the action. Instead we consider 



only several combinations which reveal the relative im- 
portance of the terms. We list in Tab. IHII the change in 
the Polaykov loop average for different combinations of 
the operators. If one sets a goal to realistically reproduce 
the value of the Polyakov loop within 1% precision then 
the A12 + B\ + C\2 effective action is the minimal ac- 
tion which is suitable. We also demonstrate the relative 
contribution of various terms in Fig. [31 where different 
two-operator effective actions are compared to the sim- 
plest one-operator action and the A\% +B% + C12 action. 
We see that the most significant improvement comes from 
the B± term, while the next significant term is Ai- 

Next we measure the couplings of the A\ 2 + B\ + Cu 
effective model at different j3 values and corresponding 
temperatures on a 40 3 x 4 lattice, see Tab. IIV1 Note 
that among the couplings measured in the confined phase 
only potential terms and the fundamental hopping term 
have significant values. Also, for comparison, we present 
in the table the couplings for the effective model with 
all 13 operators (computed at (3 — 2.7310). We use the 
same statistics for this measurement. It is obvious that 
the errors are much larger for the full action. We also 
note little change in the value of the couplings present in 
both the full and reduced A\2 + B\ + Cn effective mod- 
els. It is interesting that some of the operators present 
in the full action but not present in the reduced action 
have couplings which are significant, e.g. Xa 3 , ^b 2 1 ■ 
However this terms effect on the Polyakov loop average 
turns out to be small. 

TABLE IV: The couplings of the A12 + Bi + C12 effective 
action at various /3 and of the full 13-operator action at /3 = 
2.7310 on a 40 3 x 4 lattice. Demon measurements over 10 
sweeps after 100 sweeps of thermalization (10 runs). 



p 


1.187348 


2.2991 


2.5104 


2.7310 


2.7310 




0.0064(4) 


0.1252(3) 


0.1444(7) 


0.1666(6) 


0.1692(55) 


Abi 


-0.0436(2) 


-0.0578(2) 


-0.0679(3) 


-0.0745(3) 


-0.0742(14) 


Aa 2 


-0.0008(4) 


-0.0049(4) 


-0.0120(7) 


-0.0110(5) 


-0.0078(31) 


Aci 


-0.0002(3) 


0.0004(4) 


0.0409(4) 


0.0389(3) 


0.0225(30) 


Ac 2 


-0.0002(4) 


-0.0020(1) 


-0.0069(4) 


-0.00057(4) 


-0.0041(19) 


Aa 3 










0.0100(29) 


As 2 










-0.01377(14) 


Aa 4 










0.0013(29) 


Aa 5 










-0.0025(14) 


As 3 










-0.0037(6) 


Ac 3 










0.0008(13) 


Adi 










0.0199(29) 


Ad 2 










-0.0007(21) 



In Fig. [4]we check the action A12+B1+C12 over a wide 
range of couplings. We observe that the effective action is 
tracing the value of the Polyakov loop average computed 
from the gauge theory very closely. The agreement is very 
good in the weak coupling regime. Near the transition 
region there is small discrepancy in the values; however, 
the location of the phase transition is accurate. 
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FIG. 4: The Polyakov loop average computed from the A12 + 
Bi + C\2 effective model and the gauge theory on a 
lattice. 



IV. SUMMARY 

We presented a simple empirical method of selecting 
relevant terms in the construction of an effective spin 
action capable of reproducing a single observable (the 
Polyakov loop average). By choosing operators for the 
effective action and measuring their couplings using the 
demon method, we were able to construct and test var- 
ious effective action models. The simplest one-operator 
effective action was analyzed with respect to the size of 
non-equilibrium effects of thermalization of the expecta- 



tion value of the Polyakov loop. It was also compared to 
perturbation theory results. We find that the strong cou- 
pling regime, where the perturbation theory results cor- 
respond to the measured observable values, is the regime 
of minimal flow. The weak coupling regime exhibits seri- 
ous non-equilibrium effects, which indicates that in this 
region the one-operator action is not appropriate for the 
description of Polyakov loop configurations. More gen- 
eral effective actions with up to 13 operators were con- 
sidered and analyzed with respect to Polyakov loop flow. 
Setting as a criterion to reproduce the expectation value 
of the Polyakov loop with accuracy under 1%, we found 
that it is enough to limit the effective action to five opera- 
tors Ai2+Bi+C\2- We showed that this action is capable 
of correctly reproducing the value of the Polyakov loop 
average over a wide range of lattice couplings. 

We should note that it is also possible to consider dif- 
ferent observables. This would, in general, require a con- 
struction of a new action. It is obvious that for long 
distance observables, such as Polyakov loop correlators, 
one would need a more complex effective action model. 
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